Do baryons trace dark matter in the early universe? 
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Baryon-density perturbations of large amplitude may exist if they are compensated by dark- matter 
perturbations such that the total density is unchanged. Primordial abundances and galaxy clusters 
allow these compensated isocurvature perturbations (CIPs) to have amplitudes as large as ~ 10%. 
CIPs will modulate the power spectrum of cosmic microwave background (CMB) fluctuations — 
those due to the usual adiabatic perturbations — as a function of position on the sky. This leads to 
correlations between different spherical-harmonic coefficients of the temperature and/or polarization 
maps, and induces polarization _B-modes. Here, the magnitude of these effects is calculated and 
techniques to measure them are introduced. While a CIP of this amplitude can be probed on large 
scales with existing data, forthcoming CMB experiments should improve the sensitivity to CIPs by 
at least an order of magnitude. 

PACS numbers: 98.80.-k 



We have been conditioned to believe that the ~ 10~^ 
variations in the cosmic microwave background (CMB) 
temperature [ij, |2| imply that the matter in the early 
Universe was distributed with similarly small variations. 
This is certainly true if primordial perturbations are adia- 
batic; i.e., if there are perturbations only to the tote/ mat- 
ter content, with the fractional contributions of baryons, 
dark matter, photon, and neutrinos the same everywhere. 
It is also true for many isocurvature models [3| , where the 
total density is fixed. 

It is therefore a surprise that perturbations in the 
baryon density can be almost arbitrarily large — far larger 
than 10~^ — as long they are compensated by dark-matter 
perturbations in such a way that the total-matter density 
remains unchanged [J, |5| . These compensated isocurva- 
ture perturbations (CIPs) induce no gravitational fields, 
as the total matter density in this mode is spatially ho- 
mogeneous. Baryon-pressure gradients induce motions 
at the baryon sound speed which, at the time of pri- 
mary CMB decoupling {z ^ 1091, decoupling hereafter), 
is (v/c) - (T/mp)i/2 ^ (eV/GeV)i/2 ^ IQ-^-S. These 
motions affect the photon temperature only on distances 
< xo~^-^ times the horizon at decoupling, that is, CMB 
multipole moments with / > 10^ [J|, far larger than those 
{I < 10**) probed by CMB experiments. Thus, while 
the CMB power spectrum currently constrains the mean 
baryon-to-dark-matter ratio precisely, it tells us nothing 
about spatial variations in this ratio. 

Big-bang nucleosynthesis (BBN) and galaxy-cluster 
baryon fractions constrain the CIP amplitude to be less 
than ~ 10% [51 . Consequences of CIPs for galaxy surveys 
are small [J]. Measurements of 21-cm radiation from the 
dark ages would be sensitive to CIPs [J, la, lZ| , but such 
measurements are still off in the future. 

Here we show that the primordial relative distribution 
of baryons and dark matter can be determined with the 
CMB. Our principle motivation is curiosity — is the com- 



mon assumption that baryons trace dark matter in the 
early Universe justified empirically? However, a search 
for CIPs is also motivated by the curvaton models [8[ that 
predict their existence [J, |9[ and perhaps by recent ideas 
linking baryon and dark- matter densities [lO|. Moreover, 
if the Planck satellite finds evidence for primordial isocur- 
vature perturbations, the CIP measurements we describe 
below will be essential to determine how that perturba- 
tion is distributed between baryons and dark matter. 

CIPs modulate the baryon and dark-matter densities 
at decoupling, where ^ 90% of photons last scatter, and 
at reionization, where ~ 10% of CMB photons last scat- 
ter. There will thus be a modulation of the small-scale 
temperature and polarization power spectra from one 
patch of sky to another. This modifies the power spec- 
trum obtained by averaging over the entire sky, induces 
polarization B modes, and causes correlations between 
different spherical-harmonic coefficients of the tempera- 
ture and/or polarization maps. The effects on the CMB 
are analogous to those of gravitational lensing llj . The 
B-mode power spectrum induced by CIPs through the 
modulation of the reionization optical depth has already 
been calculated [5|. 

We show, however, that the CMB effects induced by 
modulation of the baryon density at decoupling are con- 
siderably larger than those induced at reionization. Our 
calculation follows Ref. [12| , where the CMB effects of a 
spatially- varying cosmological parameter (there the fine- 
structure constant) were considered. This variation in- 
duces a spatially-varying power spectrum. We have ex- 
tended the formalism of Ref. [13] to calculate the effect of 
CIPs on top of the usual adiabatic initial conditions, ex- 
tending the flat-sky formalism developed there to the full 
sky, and generalizing the calculation to scales of smaller 
width than the decoupling surface. Since the technical 
details are complicated, we present them elsewhere [13| 
and focus here on the principal science results. 
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FIG. 1: The B-mode power spectrum (in ^K^) for the CIP- 
induced contribution from decoupling (solid black curve), con- 
trasted with contributions from patchy scattering (short-dash 
green curve), and patchy screening (dotted red curve) at 
reionization. We use a scale-invariant spectrum of CIPs with 
the amplitude A ~ 0.017 that saturates the galaxy-cluster 
bound. 



The CIP involves baryon and cold-dark-niatter den- 
sities pbix) = pb[l + A(x)] and pc(x) ^ pc - /0&A(x), 
written as functions of position x in terms of a fractional 
baryon-density perturbation A(x). Note that the total 
matter density pfc(x)-fpc(x) associated with the CIP does 
not vary with x. We assume that A(x) is a random field 
with a scale-invariant power spectrum P^{k) = Ak^^ , as 
may be expected if CIPs arise somehow from inflation, 
and A is a dimensionless amplitude. The rms variation 
Aci in the baryon-to-dark-matter ratio between galaxy 
clusters obeys the constraint Ad ^ 0.08. 

When the three-dimensional field is projected onto a 
narrow spherical surface, the resulting angular power 
spectrum for A will be C^ ~ A/(7ri^) for mulipole mo- 
ments L <, (ryo — rj\s)/ (Jr,, where ryis and 7]q arc the confor- 
mal time at last scatter and today, respectively, and cr,, is 
the rms conformal-time width of the last-scattering sur- 
face. At smaller angular scales (larger i), the variation 
in A is suppressed by the finite width of the scattering 
surface. The angular power spectrum for A can then 
be approximated by C^ ~ ^(^o — VIb)/{'^\/^L^<^t]) for 
L <: (770 — Vis) /cTri [si ll3|. The rms variation Ad in the 
baryon-to-dark-matter ratio on galaxy cluster scales is 

A?i = ^ I fc2 dk [3n{kR)/ikR)f PA(fc), (1) 

where R is the mean separation between galaxy clus- 
ters. The integral has a formal logarithmic divergence 



at low k which is cut off, however, by the horizon 
fcmin ~ (10 Gpc)"^ Taking i? ~ 10 Mpc, we find A^j ~ 
ylln(1000)/27r2. Thus, Ad < 0.08 implies A < 0.017. A 
weaker bound {A < 0.046) comes from BBN. 

Now consider the CMB fluctuations produced at de- 
coupling. The CIP-induced variation of the baryon and 
dark-matter densities across the sky modulates the small- 
scale power spectrum, and this modulation induces off- 
diagonal correlations in the CMB [12|, ll3| . 

Moreover, B modes are induced in the CMB polar- 
ization |12| . The induced spherical-harmonic coefficients 
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where A^m arc the spherical-harmonic coefficients for 
A(n); daj^^,/dA is the derivative of the usual _E-mode 
spherical-harmonic coefficient with respect to A (com- 
puted using the CAMB code [l^) and 
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This induced B mode arises by modulating the first-order 
adiabatic perturbation to first order in the CIP. This is 
because the sound speed, photon diffusion length, and 
visibility function, assumed spatially homogeneous in the 
standard treatment, all depend on the local baryon den- 
sity. In contrast, when the CMB is gravitationally lensed, 
daf^, /dA is replaced by a function encoding a deflection. 

Figure [T] shows the results of our calculations for 
the i3-mode power spectrum C^^ induced by a scale- 
invariant spectrum of CIPs with the largest amplitude 
{A ~ 0.017) consistent with galaxy-cluster baryon-to- 
dark-matter ratios. CIPs modulate the reionization opti- 
cal depth, and as noted in Ref. [Sj, this also generates B 
modes, through patchy screening [1J| and scattering [15[ 
of primordial CMB fluctuations. We plot these reion- 
ization contributions in Fig. [T] for the same CIP am- 
plitude. We see that the i?-mode power spectrum in- 
duced at decoupling is larger (by up to 3 orders of mag- 
nitude) than that induced at reionization, for I ^ 50. The 
decoupling-induced B modes are larger because (a) they 
involve ^ 90% of the photons, rather than ^ 10%, and 
(b) the flnite width of the reionization re-scattering sur- 
face smooths the angular A fluctuations to larger angular 
scales (lower L) than it does for decoupling. Reconstruc- 
tion of Alm depends primarily on higher-/ modes, and 
so the baryon-density modulation at decoupling is more 
important in probing CIPs than that at reionization. 

Wc now turn to the reconstruction of CIPs from 
CMB maps. In the absence of CIPs, the multipole mo- 
ments af^ obey the relation (af^a^'^.j = C^^'SwSmm' 
{X G {T,E,B}); i.e., spherical-harmonic coefficients 
with (Im) 7^ {I'm') are statistically independent. How- 
ever, if there is spatial modulation of the power spectrum. 
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FIG. 2: Shown are the errors in C^ from the TT, TE, EE, TB, and EB estimators for the CIP perturbation A at the surface 
of last scatter for (a) WMAP and (b) a CMB polarization satellite, with the specifications spelled out in the EPIC mission 
concept study. Also shown (signal) is the power spectrum C^ for a scale-invariant spectrum of CIPs with the maximum 
amplitude allowed by galaxy clusters. 
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FIG. 3: (a) Shown are the total expected errors in C^ from the combined TT ^ TE, EE, TB, and EB estimators for the CIP 
perturbation A at the surface of last scatter for several current and forthcoming CMB experiments. Also shown (signal) is 
the power spectrum C^ for a scale-invariant spectrum of CIPs with the maximum amplitude allowed by galaxy clusters, (b) 
The signal-to-noise as a function of the rms fluctuation in the galaxy cluster baryon-to-dark-matter ratio along different lines 
of sight at 2 ~ 1091. The vertical line shows the upper limit to the rms amplitude from galaxy clusters (excluded region is 
shaded) . 



then there will be ofF-diagonal (/ ^ I' , m ^ m!) correla- 
tions, 

+ Y^D^M,xx'^^M^^^^ (4) 
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where D,,, ' = Ai^mS,,,' are bipolar spherical 
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harmonics [16[ , and 5*;;,' are coupling coefficients that 



are calculated in Ref. [13|. 

Construction of minimum- variance estimators for the 



Alm and their associated errors is straightforward [17 
For example, for the EB correlation, the estimator is 
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where C^-^ are power spectra including noise, and 
Djj, IS the minimum-variance estimator tor 1^;^, 

|17| . From these one can estimate C^ . 

Figure [2] shows the predicted errors in the CIP power 
spectrum reconstruction from the TT, EE, TE, TB, and 
EB estimators for the Wilkinson Microwave Anisotropy 
Probe (WMAP) and the proposed Experimental Probe 
of Inflationary Cosmology (EPIC). These are SCf-^ = 
[i2L + 1)]-^/' (^if )V/sky, where the (2L + l)-i/2 fac- 
tor results from the multiple modes available at each L. 

Instrumental parameters for WMAP are a beamwidth 
of 21' (full- width half-max), noise-equivalent tempera- 
ture (NET) of 1200 fJ.Ky/s, fraction of sky analyzed 
/sky = 0.65, and observation time iobs = 7 years. For 
EPIC (150 Ghz channel) we assume a beamwidth of 5 
arcmin, NET of 2.0 /iK-^s, and observation time fobs = 
4 years, also with /sky = 0.65. 

For WMAP, the best sensitivity comes from TT. For 
EPIC, the sensitivity at L > 100 comes primarily from 
the TB estimator. We have checked that the best sen- 
sitivity for Planck comes from TT, while some ground- 
based experiments (e.g., SPTPol) benefit from polariza- 
tion. 

The left panel of Fig. [3] shows the errors in the CIP 
power spectrum reconstruction obtained by combining 
the TT, TE, EE, TB, and EB estimators for a variety 
of CMB experiments. The signal to noise ratio is given by 
S/N = {ifs,.y/2)j: u2i2L+l)[iCt/cTlJ]'y/\ 

^ J sky 

The right panel shows the S/N for detection of a scale- 
invariant spectrum of CIPs as a function of Ad. With 
WMAP, a CIP saturating the cluster bound is marginally 
accessible on the largest scales. Planck should be able to 
probe rms CIP amplitudes of 3 x 10^^ and higher. Sig- 
nificant improvements in sensitivity should be obtained 
with upcoming experiments like Polarbear, SPTPol, and 
ACTPol. We see that S/N values > 3 may be possible 
with EPIC for an rms CIP amplitude of 4 x 10~^, a factor 
of ~ 20 lower than the current limit. 

The tools for these measurements should be gener- 
alizations of those used for weak Icnsing of the CMB 
[ISj . which also produces off-diagonal correlations. CIPs 
should be distinguishable from lensing, since these phys- 
ical effects are distinct, as evidenced by differing forms 

T V V' 1 L 

for the coupling coefficients S^l . In Ref. [19[, it is 
shown that for the analogous case of patchy reionization, 
optical depth fluctuations may be separated from lensing 
without significant loss in S/N , and we expect that this 
is also true for CIPs. We leave for future work the de- 
velopment of tools to distinguish CIPs from weak lensing 
and contaminants like Galactic foregrounds. 



CIPs are an intriguing possibility and a prediction of 
some inflationary models. With the measurements we 
have described here, we may soon know empirically how 
closely dark matter and baryons trace each other in the 
early Universe. 



We thank G. Holder, T. L. Smith, M. Lo Verde, C. Chi- 
ang, K. M. Smith, M. Zaldarriaga, and D. N. Spergel for 
stimulating conversations. We thank B. Jones and A. 
Fraisse for providing Spider parameters. DC was sup- 
ported by NSF AST-0807044 and MK by DoE DE-FG03- 
92-ER40701 and NASA NNX10AD04G. Part of this re- 
search was carried out at the Jet Propulsion Laboratory, 
Caltech, under a contract with NASA. 



C. L. Bennett, A. Ban day, K. M. Gorski e t al., Astrophys. 
J. 464, L1-L4 (1996) ■astro-ph/9601067' . 
E. Komatsu et al, Astrophys. J. Suppl. 192, 18 (2011) 
[arXiv: 1001. 4 538 . 

R. Bean, J. Dunkley and E. Pierpaoli, Phys. Rev. D74, 
063503 (2006) astro-ph/0606685 . 

C. Gor don and J. R. Prit chard, Phys. Rev. D 80, 063535 
(2009) | arXiv:Q907.5"400] . 

G. P. Holder, K. M. NoUett and A. van Engelen, Astro- 
phys. J. 716, 907 (2010) arXiv:0907.3919;. 
R. Barkana and A. Loeb, Mon. Not. Roy. Astron. Soc. 
Lett. 363, L36 (2005) ast ro-ph/050 2083 . 
A. Le wis and A . Challinor, Phys. Rev. D 76, 083005 
(2007) 'astro-ph /0702600| . 



D. H. Lyth, C. U ngarelli and D. Wa nds, Phys. Rev. D 

67, 023503 (2003) astro-ph/0208055 . 
[9] C. Gordon and A. L ewis, Phys. Rev. D 67, 123513 (2003) 

[astro-ph/0212248| . 
[10] D. E. Kaplan, M. A. Luty and K. M. Zurek, Phys. Rev. 

D79, 115016 (2009) [arXiv: 090 1.4117] . 
[11] A. Lewis and A. Challinor, Phys. Rept. 429, 1 (2006) 

astro-ph/0601594 . 
[12] K. Sigurdson, A. Kurylov and M. Kaniionkowski, Phys. 

Rev. D 68, 103509 (2003) astro-ph /0306372,. 

[13] D. Grin, O. Dore, M. Kamionkowski. [arXiv:1107.5047| . 
[14] C. Dvorkin an d K. M. Smith , Phys. Rev. D 79, 

043003 (2009) [arXiv:0812.1566] : C. Dvorkin, W. Hu 

and K. M. Smith, Phys. Rev. D 79, 107302 (2009) 
arXiv:0902.4413 . 
[15] W. Hu, Astrophys. J. 529, 12 (2000) [astro-ph/9907103] ; 

D. Baumann, A. Cooray and M. Kamionkowski, New 

Astron. 8, 565 (2003) [astro-ph/02085TT ; . 
[16] A. Hajian and T. Souradeep, Astrophys. J. 597, L5 

(2003) astro-ph/0308001 . 
[17] A. R. PuUen and M. Kamionkowski, Phys. Rev. D 

76, 103529 (2007) arXiv:0709.1144 ; M. Kamionkowski , 

Phys. Rev. Lett. 102, 111302 (2009) [arXiv:0810.1286] : 

V. Gluscevic, M. Kamionkowski and A. Cooray, Phys. 

Rev. D 80, 023510 (2009) arXi v:0905.l687] . 
[18] K. M. Smith, O. Zahn and O. Dore, Phys. Rev. D76, 

043510 (2007) arXiv:0705 .3980 ; C. M. Hirata et al, 

Phys. Rev. D78, 043520 (2008) [arXiv: 080 1.0644] : S. Das 

et al, arXiv:1103.2124.. 



[19] M. Su, A. P. S. Yadav, M. McQuinn, J. Yoo, M. Zaldar- 
riaga, .arXiv:1106.4313j . 



